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Abstract 

o 
o 

1-1 Motivated by the recent achievements in manipulation of C6o molecules in STM exper- 

iments, we study theoretically the structure and electronic properties of a C6o molecule in 
an STM-tunnelj unction with a magnetic tip and magnetic adatom on a Cu(lll) surface from 

OO 

first-principle calculations. For the case of V tip/adatom, we demonstrate how spin-coupling 
between the magnetic V atoms mediated by the C6o can be observed in the electronic trans- 
port, which display a strong spin-filtering effect, allowing mainly majority-spin electrons to 
pass(>95%). Moreover, we find a significant change in the conductance between parallel and 

X 

anti-parallel spin polarizations in the junction (86%) which suggests that STM experiments 
should be able to characterize the magnetism and spin-coupling for these systems. 
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INTRODUCTION 



Organic materials typically offer small spin-orbit and hyperfine interactions, which are prerequi- 
sites for spintronic applications, due to allowing long spin lifetimes. Thus there is a great interest in 
organic building blocks for molecular spintronics J^l an d a thorough understanding of spin trans- 
port and magnetism in these systems is called for. It is therefore important to establish model 
molecular spintronic systems where spin-transport and magnetic interactions can be examined ex- 
perimentally. 

Recently, it has been demonstrated in low temperature scanning tunneling microscopy (STM) 
experiments how C6o molecules can be picked up by the STM-tip, and subsequently controllably 
used to contact structures such as adatoms, clusters, and molecules placed on a substrate surface.'^' 
The Ceo molecule is considered as an attractive anchoring group for molecular electronics due to 
its mechanical robustness. 8 Moreover, the lowest unoccupied molecular orbital (LUMO) of C60 is 
close to the Fermi level of ferromagnetic elements, which makes spin injection relatively easy,^ 
and makes C6o a promising building block in molecular spintronics. 

The high symmetry of the C6o allows detailed characterization of the bonding geometries in 
STM. In particular, one can determine which part of C6o is pointing towards the tip/surface prior 
to contact formation, and also afterward, while the C(,o is placed on the tip. 7 Subsequently, it is 
possible to investigate the interactions between tip and sample via electronic transport measure- 
ments, as the tip and and sample are brought into contact. STM is a powerful tool for investigating 
spin-transport in magnetic nanostructures.^^ Direct magnetic interaction between STM tip and 
magnetic materials on a substrate has been studied in number of works ,E2Hm and STM has been 
used to probe spin in organic molecules.^In the case of a magnetic tip and magnetic surfaces, this 
method may be used to study spin-transport and interactions though organic molecular systems 
bound to the surface and gain insight into single-molecular magnetic properties. 

Among organic compounds, carbonic rings when combined with transition metals show versa- 
tile combinations for molecular spintronics.^ The interaction of n electron system of such rings 
with d orbitals of the transition metals, is a key to electron and subsequently spin transport. One 
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example of such systems is presented in a theoretical study, where calculations have been used 
to examine spin transport in a benzene-Co system on a Cu(OOl) surface contacted by a Cr tip.^1 
Further on, the magnetic properties and spin transport have been calculated for organometallic 
"multidecker" wires where magnetic atoms are sandwiched between organic parts. 21 Multidecker 
systems involving vanadium appear very promising due to their half-metallic behavior resulting 
in high spin-polarization of the transport.^^ Interstingly, due to the different symmetries of the 
C6o , it is possible to vary the electronic and magnetic properties depending on whether pentagon, 
hexagon or edge sites of C6o is in contact with magnetic ligand atoms. 

Here, we employ first principles calculations to predict the spin transport through spintronic 
model system consisting of a C6o molecule contacted by magnetic atoms in a STM setup. In 
particular, we predict that vanadium is a magnetic material where pronounces spin-filtering and 
spin- valve effects in STM experiments. 

SYSTEM AND METHOD 

In order to mimic a concrete STM experiment, we investigate the specific setup shown in [TJ The 
bulk regions of the contacts (i.e. STM tip and substrate) have been chosen to be non-magnetic 
copper, which has previously been employed in manipulation experiments.'^' We imagine that 
magnetic atoms are deposited on the Cu surface prior to deposition^ of the C6o molecules, and 
that the tip-electrode is prepared prior to the contact by either indenting a Cu tip into a cluster of 
these atoms in order to cover its outermost part with these, or by creating the tip from from the bulk 
magnetic material. We model the outermost part of the tip by a pyramid-like structure consisting 
of four atoms on one electrode, simulating the magnetic STM tip which is used to pick up C6o for 
subsequent contact formation to an isolated adatom on the copper (111) substrate. In the following 
we will show that the case of vanadium is remarkable. 

We used spin-polarized pseudopotential DFT calculations with the SIESTA code. 25 Electronic 
structures have been calculated within a GGA-PBE approximation to the exchange and correlation 
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functional.^ Double- £ polarized basis sets with grid cutoff of 250 Ry have been used. Spin po- 
larized transport was subsequently calculated using the non-equilibrium Green's function (NEGF) 
formalism 27 in the limit of zero voltage. In order to eliminate basis set superposition error (BSSE) 
present in methods with atomic orbitals basis sets, we have used plane-wave (PW) formalism as 
implemented in Ref. 28 for total energy calculations. In these calculations, we have been using 
ultrasoft pseudopotentials, with 30/300 Ry cutoff for wavefunction/charge density. 




Figure 1 : STM system picking up C6o with a magnetic tip, approaching a magnetic adatom on 
the nonmagnetic copper surface. The orange, red and gray spheres depict copper, vanadium and 
carbon atoms respectively. 

To guide us in our search for interesting magnetic metals to contact the C6o , and to select 
relevant bonding sites, we first performed calculations of the binding energy of a single adatom 
from the first row transition metals (from vanadium to copper) with C^q . Here our focus, besides 
achieving high magnetization, is to obtain a high binding energy in order to achieve a stable contact. 
The results are summarized in [T] for different sites on the C6o molecule. One observes that nickel 
has the strongest binding energy, but with zero total magnetization (My), and thus is probably not 
interesting for investigations of spin-transport. On the other hand, chromium enjoys the largest 
Mj, due to its largest unpaired electronic configuration [Ar]3d 5 4s 1 , but with the least binding 
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energy strength. It has already been shown that copper STM tips can pick up Noting that 

the maximum of binding energy for copper to C6o is ~ 0.8 eV, we conclude that the same action 
should be possible with V with a decent Mt. 

Table 1: Binding energy E b = E C60+a datom - E c 60 E at j a t om and total magnetization Mt — f v (n up 
n down) dr 3 per unit cell. The distance between the adatom and C6o has been optimized. 



adatom 56 (eV) 66 (eV) hexagon (eV) pentagon (eV) top (eV) Mt 



V 


-1.13 


-1.25 


-0.93 


-1.40 


-1.06 


5.00 


Cr 


-0.58 


-0.46 


1.77 


2.74 


-0.49 


6.00 


Mn 


-0.39 


-0.58 


-0.25 


-0.33 


-0.40 


5.00 


Fe 


-0.83 


0.08 


-0.87 


-0.55 


-0.75 


4.00 


Co 


-0.90 


-1.22 


-0.87 


-0.88 


-0.80 


1.00 


Ni 


-1.48 


-1.64 


-1.24 


-1.44 


-1.26 


0.00 


Cu 


-0.79 


-0.60 


-0.08 


-0.34 


-0.74 


0.00 



66 refers to sites above edges shared between two hexagons, 56 referes to sites above edges shared between a 
hexagon and a pentagon and top means top of a carbon atom. 

Based on the data of the simple guiding calculations we subsequently chose vanadium in the 
full C6o -contact simulations on the system depicted in [T] In the following we demonstrate that 
this choice for contact material in an STM setup indeed is successful in achieving a good spin- 
filter and spin-valve performance. Here, ferromagnetic (FM) and anti-ferromagnetic (AFM) spin 
alignment between atoms of tip and adatom, have been considered. The site on the C(,q with highest 
binding energy for V adatom is T] 5 , roughly over the center of a pentagon of a Cgo , and due to 
the symmetric structure of the C6o , this site is contacted both by the tip and the adatom. We find 
a binding energy of the C6o to the V-tip (upper part in [T]) is 1.3 eV, while binding of C6o to the 
adatom on the Cu substrate (lower part in [T]) is 1.1 eV. The spin-resolved transmissions for the 
FM and AFM cases are shown in|2j We first focus on the highly conducting contact configuration 
where the atomic structure of the C6o along with vanadium atoms and first copper layers of both 
sides have been relaxed to the force threshold of 0.05 eV/A. We also show the transmission spin 
polarization (TSP), defined as 



and channel decomposed transmission values in [2j Here, we first point out a remarkable spin- 
filtering effect in the FM arrangement, whereby two almost open channels conduct in the vicinity 
of the Fermi level for the majority spin component, while the minority channels are almost closed. 
For the minority spin component the resonance peaks at ~ 0.2 and 0.4 eV produce dips in the 
corresponding TSP curve, however, these will only be of importance for a voltage bias comparable 
to these energies. Transmission eigenvalues of first three dominant channels are shown in the 
third and forth panels, that clearly show two distinct channels for the FM-majority spin channels. 
Furthermore, it is striking that the channels in the AFM configuration are almost closed, except for 
small resonance peaks at ~ 0.35 eV, which again only will come into play for higher voltages. 

To better understand the nature of spin transport in the system, we have furthermore calculated 
the spatially-resolved scattering states in the contact region.^ The results are shown in [3} Here, 
we consider the conducting FM arrangement, and focus on the two eigenchannel scattering states 
with highest transmission at ^(moving in the direction up-to-down), which both are almost fully 
transmitting. For the majority spins, we notice the and d yz orbital nature of wavefunctions on 
the V adatoms contacting C(,q (z chosen perpendicular to the surface). This is in accordance with 
the Mulliken population analysis of the majority spin states of the V tip and adatoms, where the 
d u , d yz each appears half-filled. On the other hand, the J^, d x i_ y i and d\ are closer to being filled, 
while the s is closer to being empty. This points to a charge transfer from the V atoms to the C60 
leaving the d^/dyz orbital energies closest to Ep. Since the d^/dyz orbitals match the symmetry 
with angular momentum m = 1 for rotation around the V-Cgo -V axis of the pentagon-prone 3-fold 
degenerate LUMO states {t\ u symmetry! 30 * 31 } ), we can expect the observed orbitals in the transport 
channels. For minority spins d^, d yz orbitals are almost empty and shifted away from Ep, resulting 
in a vanishing transmission. The rotational symmetric m = channels appear as resonances in the 
channel transmissions above Ep and thus play a minor role. 

In typical STM experiments the conductance is probed from the tunnel-regime to contact. We 
have performed transport calculations as the tip is approaching the surface adatom, until the tip- 
molecule distance (d shown in[T]) approximately reaches the equilibrium distance discussed above. 
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Figure 2: Transmission spectra for FM and AFM arrangements. The first row shows spin-resolved 
transmission spectra for each arrangement. The second row demonstrates the corresponding trans- 
mission spin polarization as defined in the text. The third and forth rows show transmission eigen- 
values for three channels in each spin manifold. 
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Figure 3: Scattering states at E = Ep of first two dominant eigenchannels for (a,b) majority and 
(c,d) minority spin components in FM arrangement. Blue and red indicate the positive and negative 
sign of the real part of the wavefunction. 
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In [4] we display the conductance along with the corresponding TSP as a function of tip distance. 
One can observe the trend of increasing the conductance of the majority spins and thus TSP, in 
the FM case, while in AFM, the conductance values are considerably smaller all the way to the 
equilibrium contact distance. 



FM AFM 




Figure 4: Spin-resolved conductance and transmission spin polarization(TSP) vs. C6o -adatom 
separation. 

The marked difference between the FM and AFM conductance properties indicates how it is 
possible to probe the spin-coupling mediated by the Cgo between the magnetic tip and substrate. 
The calculated magnetic interaction between the tip and adatom, the magnetic exchange energy 
defined as Epm — Eafm, is shown in^a) when the tip-molecule is approaching the adatom. This 
shows that the FM arrangement becomes favorable as the molecule reaches the equilibrium dis- 
tance to the surface adatom. To be sure about the fidelity of the values obtained here, we have 
performed the same study using the PW method. We find that the trend is the same and the values 
are even more pronounced in favor of FM arrangement, though of the same order of magnitude. 

A graph of the total conductance versus C6o -adatom distance is shown in [5jb). The con- 
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ductance behavior demonstrates a magnetic valve, being closed for FM and open for AFM, if 
we imagine an external control over the magnetization of tip/substrate. In a typical experiment 
with a bulk magnetic tip the magnetization of the tip will be determined by the intrinsic magnetic 
anisotropy of the crystalline magnetization, which fix the magnetization axes. As the tip-molecule 
approaches the adatom on the non-magnetic surface, the magnetization of this will be determined 
by the interaction with the tip mediated by the molecule. In this case the adatom magnetization 
will align according to the thermal occupations. Using the calculated energies for the FM/AFM 
configurations we can now calculate the conductance for a given temperature for each distance, 

e- E ™( d V k * T G FM (d) + e- E ^ d V k * T G AFM (d) 

[ ' > e -E FM (d)/k B T +e -E AFM (d)/k B T ■ (Z) 

This is shown in^c) for temperatures of 10K, 50K and 300K, showing only the FM configuration 
at low temperature. 

The absolute distance is typically not known in an actual STM experiment. In principle, a par- 
ticular conductance could be realized with both FM or AFM spin configurations - a conductance 
of e 2 /h could result from a single spin-channel with perfect transmission or two half-transmitting 
channels. In combination with measurements of the conductance, measurements of current shot- 
noise as characterized by the Fano factor, 

F _ L,ME F )(l-T n (E F )) ^ 

Y. n Tn(E F ) 

can provide further insights into the distribution of transmissions in the conductance channels, as 
demonstrated for molecular contacts j n ^d), we observe how the noise is significantly smaller 
for the FM configuration and drops already well before contact (d — d eq ~ 1.5 A) is established. It 
can also be inferred how the transport is carried by two channels in the case of FM in contact. 

In summary, we have shown how C$q in a STM junction formed between a V tip and a V 
adatom on a Cu substrate display highly spin-polarized transport and a remarkable difference in 
conductance and noise properties depending on the spin- alignment of the tip and adatom. Our 
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Figure 5: (a) magnetic exchange energy, (b) conductance for FM and AFM configurations, (c) 
thermal average of the conductance for different temperatures, (d) Fano factor of transmission as a 
function of C6o - adatom separation for the FM and AFM configurations. 
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results illustrate how the unique manipulation capabilities of the STM and the high symmetry of 
the C6o may be exploited in order to gain insights into spin-transport and magnetic interactions at 
the molecular scale. 
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